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ABSTRACT 

Observations of AGN have suggested a possible delay between the peak of star formation (on 
some scale) and AGN activity. Inefficient fueling (and/or feedback) from fast stellar winds 
has been invoked to explain this, but we argue this is unlikely in bright systems accreting 
primarily cold dense gas. We show that such a delay can arise even in bright quasars for 
purely dynamical reasons. If some large-scale process produces rapid inflow, smaller scales 
will quickly become gas-dominated. As the gas density peaks, so does the SFR. However, 
gravitational torques which govern further inflow are relatively inefficient in gas-dominated 
systems; as more gas is turned into stars, the stars provide an efficient angular momentum 
sink allowing more rapid inflow. Moreover, the gas provided to the central regions in mergers 
or strong disk instabilities will typically be ^ 100 times larger than that needed to fuel the 
BH; the system is effectively in the "infinite gas supply" limit. BH growth can therefore 
continue for some time while the gas supply exhausts, until it has significantly depleted to 
the point where the BH is finally "starved." Both of these effects act together with comparable 
magnitude, and mean that the peak of BH growth can lag the peak in the SFR measured at a 
given scale by a timescale corresponding to the gas exhaustion time on that scale 10—100 
local dynamical times). This predicts that the inferred delay will vary in a specific manner with 
the radius over which the star formation rate is measured. We discuss possible implications 
for the role of AGN feedback in suppressing star formation activity. 

Key words: galaxies: formation — galaxies: evolution — galaxies: active — quasars: general 
— cosmology: theory 



1 INTRODUCTION 



The idea of a connection between AGN and starburst activity 
in galaxies has a long history, going back to speculation that 
dense nuclear gas concentrations in ultra-luminous infrared galax- 
ies (ULIRGs) fuel super-massive black holes that eventually appear 
as luminous quasars, expelling the nuclear gas and dust jSanders| 
|1999[[Jose ph 1999'). The discovery of tight correlations between the 
mass of nuclear black holes and their host bulge properties (IMagor-| 
naiTet al. 1998, Ferrarese & Merritt 2000; Gebhardt et al. 2000 



Hopkins et al.,,2007] [Aller & Richstone||2007t feoh & Mancini| 



2009 1 has suggested that black hole and galaxy formation are co- 



eval, and the cosmic history of black hole growth has roughly sim- 
ilar shape to the cosmic star formation history (see e.g. jMerlonij 
|et al.||2004| >. Observations have found AGN to preferentially live 
in hosts with enhanced star formation activity relative to control 
galaxy samples (|Brotherton et al. 1999, Canalizo & Stockton|2001J 
Kauffmann et al. 2003, Jahnke et al. 2004, Sanchez et al. 2004 
Vanden Berk et al.,2006, Kaviraj et al.,201 1, Zakamska et al.,2006^ . 



However, some recent observations have suggested that AGN 
and starburst activity may not be contemporaneous, even in a sta- 
tistical sense. Specifically, |Davies et al.|(f2007i study a local sam- 
ple of AGN at high resolution and see tentative suggestions that 
strong AGN are not present until the stellar populations in the cen- 
tral < 10— lOOpc are of order a few lO'yr old. On large scales, 
Schawinski et al. ( 2009') and ' Wild et"aL] ( [20T0l l (see also |Trichas| 
et al..2009, ,Tadhunter et al.,20TT] ) compare the level of AGN ac- 
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tivity to the integrated (galaxy-wide or central-kpc, respectively) 
stellar populations of local galaxies, and suggest that AGN activity 
may increase after the global stellar populations reach an age of a 
few 10* yr 

[Norman & Scovillel ( |1988| l suggested that a delay between SF 
and AGN activity might result from stellar evolution timescales if 
the BH is fueled directly by stellar winds. At low black hole accre- 
tion rates (BHARs), this may be important: young stellar popula- 
tions (< lOOMyr) produce fast but tenuous O-star winds, which 
contribute to a diffuse ISM that is hot and/or has large turbu- 
lent/bulk motions of ^ 100km s~ ' . If the BHAR is dominated by 
the accretion of gas from this hot/diffuse phase of the ISM, then 
the Bondi-Hoyle accretion rate is suppressed by a factor oc c~ff , 
where Ccff is the effective sound speed and/or turbulent speed in the 
ISM. As a consequence, if this energy cannot be efficiently radi- 
ated, larger accretion rates will require the arrival of a higher frac- 
tion of "slow winds" from more evolved stellar populations, which 
(having lower velocities and higher densities) can be more easily 
accreted. A stronger extension to this argument might even be that 
the energy/momentum feedback in these winds suppresses accre- 
tion even of colder gas for some additional time (Wild et al. 201^. 
There are suggestions that accretion at very low Eddington ratios, 
may be dominated by accretion of diffuse/hot gas, in 
which case the nature of stellar ejecta could have dramatic effects 
l |Allen et al.|2006l|Soria et al.|2006||Hrckox et al.|2009| ). 

However, at high accretion rates > 1% of Eddington, it is gen- 
erally believed that BHs are accreting from a dynamically cold, 
disky high-density medium. For example, maser observations have 
shown typical material at ~ 1 — lOpc has densities > lO^cm"^, 
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with effective te mperatures < lOOK and turbulent velocities of 
~ 10 -50km s"' (Greenhill et al. 2003 : Lonsdale et al.'2003"Fr^ 
uscione et ari|2005t Henkel et al. 2005 ; Kondratko et al. 2006. 



Nenkova et al.| 2008 ^ . The cold, disky, dense phase ex- 
ists in e.g. dense starburst systems on ~ 100— lOOOpc scales re- 
gardless of the presence of extremely young stellar populations 



iDownes & Solomon|1998i|Bryant & Scoville|1999||Tacconi et al.| 



1 1999 1. This is not surprising: tenuous winds will have little dynam- 
ical/accelerative effect on cold clouds at these densities, and the gas 
is believed to be driven in from larger radii, not locally produced 
by stellar ejecta. For typical conditions, with M, ~ Mgas on scales 
< lOpc, the mass/momentum injection rates from O-star winds are 
much smaller than the inflow rate/gravitational force (for a nuclear 
SFR similar to observed values in IDavies et al.ll2007l and wind 
momentum-loading from |Leitherer et al.|1999^ . Moreover, even if 
the cold medium is replenished or heated by winds, the cooling 
time in such regions is ~ 10~' times shorter than the dynamical 
time, so the medium can immediately re-form a cold disk (it retains 
little memory of the wind properties; see e.g. [Schartmann et al.] 
[20T0l >. 

In this paper, we therefore investigate whether "delays" simi- 
lar to those observed can appear independent of stellar fueling (or 
feedback) processes. We show, using a series of high-resolution hy- 
drodynamic simulations, that basic dynamical processes can give 
rise to qualitatively similar delays even in AGN accreting from 
cold, dense media. 

2 THE SIMULATIONS 



[Hopkins & Quataert|p010| l give a detailed description of the simu- 
lations used here; we briefly summarize some important properties. 
The simulations were performed with the parallel TreeSPH code 
GADGET-3 jSpringel|2005^ . They include stellar disks and bulges, 
dark matter halos, gas, and BHs. 

Because of the large dynamic range in both space and time 
needed for the self-consistent simulation of galactic inflows and nu- 
clear disk formation, we use a "zoom-in" re-simulation approach. 
This begins with a large suite of simulations of galaxy-galaxy 
mergers, and isolated bar-(un)stable disks. These simulations have 
0.5 X lO' particles, corresponding to a spatial resolution of 50pc. 
These simulations have been described in a series of previous pa- 
pers (IDi Matteo et al."2005''Robertson et a l.|2006b||Cox et al.|2006| 
[hung er et al. 2008, Hopkins et al. 2009 ). 

Following gas down to the BH accretion disk requires much 
higher spatial resolution than is present in the galaxy-scale simu- 
lations. We therefore select snapshots from the galaxy-scale sim- 
ulations at key epochs and isolate the central 0.1 — 1 kpc region, 
which contains most of the gas that has been driven in from large 
scales.' From this mass distribution, we then re-populate the gas in 
the central regions at much higher resolution, and simulate the dy- 
namics for many local dynamical times. Specifically, we either take 
the mass distribution in e.g. the central kpc "as is" from the par- 
ent distribution ("cutting out" the central region and re-populating 
it with higher-resolution particles), or fit it to an exponential disk 
model which we initialize as an idealized BFl/disk/bulge system 
in equilibrium (to avoid artificial features from the low-resolution 
initial simulation). These simulations involve 10* particles, with a 
resolution of a few pc and particle masses of 10"* Mq. We have 
run ~ 50 such re-simulations, corresponding to variations in the 
global system properties, the model of star formation and feedback, 



Typically ~ 10"'A/q of gas, within a scale length of ~ 0.3 — 0.5kpc. 



and the exact time in the larger-scale dynamics at which the re- 
simulation occurs. [Hopkins & Quataert| ( |2010 1| present a number of 
tests of this re-simulation approach and show that it is reasonably 
robust for this problem. This is largely because, for gas-rich disky 
systems, the central ~ 300 pc becomes strongly self-gravitating, 
generating instabilities that dominate the subsequent dynamics. 

These initial re-simulations capture the dynamics down to 
~ 10 pc, still insufficient to quantitatively describe accretion onto 
a central BH. We thus repeat our re-simulation process once more, 
using the central ~ 10 — 30pc of the first re-simulations to initialize 
a new set of even smaller-scale simulations. These typically have 
~ 10* particles, a spatial resolution of 0.1 pc, and a particle mass 
~ 100 Mq. We caiTied out ^ 50 such simulations to test the robust- 
ness of our conclusions and survey the parameter space of galaxy 
properties. These final re-simulations are evolved for ~ 10^ years. 
We also carry out a few extremely high-resolution intermediate- 
scale simulations, which include ~ 5 x 10^ particles and resolve 
structure from ~ kpc to ~ 0.3 pc, obviating the need for a second 
zoom-in iteration. 

Because of the one-way "re-simulation" method adopted here, 
the small-scale simulations are effectively instantaneous relative to 
the larger scale simulations. They capture the gas dynamics for 
one particular realization of the gas distribution on small and large 
scales, but cannot be evolved for long timescales relative to the 
large-scale simulations. 

Our simulations include gas cooling and star formation, with 
gas fo rming stars at a rate motivated by the observed IKennicutt] 
( 1998| relation (p, oc p''^). Because we cannot resolve the detailed 
processes of supemovae explosions, stellar winds, and radiative 
feedback, feedback from stars is modeled with an effective equation 
of state ( [Springel & Hemquist|200 3 ). In this model, feedback is as- 
sumed to generate a non-thermal (turbulent, in reality) sound speed; 
we use sub-grid sound speeds ~ 20 — lOOkms"', motivated by a 
variety of observations of dense, star forming galaxies (|Downes &. 
Solomon 1998; Brya nt & Scoville|[T999i [Forster Sch reiber et"d! 
|2006t |Iono et al.,2007) . Within this range, we found little differ- 
ence in the physics of angular momentum transport or in the result- 
ing BHARs, gas masses, etc. ( [Hopkins & Quataert|2010[ l (see also 

Critically, our model does not include any explicit "fast stellar 
winds" or time-dependence in the stellar feedback or mass recy- 
cling; there is an effective ISM pressure which corresponds to tur- 
bulent ISM sound speeds < 100km s~', but this is much less than 
the ^ 1000 — 2000 km s~ ' O-star wind velocities and, more impor- 
tantly, is time-independent. Therefore any characteristic timescale 
that emerges in these simulations is not a consequence of stellar 
evolution or mass recycling (as we have implemented it), even if 
this is in reality critical to explaining certain properties of the ISM 
(see e.g. [Hopkins et al.|201 la c bl. 

For all simulations, the total SFR within a given radius is di- 
rectly determined by integration of the SFR in each gas particle 
(this is well-converged in each case with respect to resolution; see 
[Springe l et aH2005|l. The BHAR requires more care. In the nuclear- 
scale simulations, the BHAR (which we assume is proportional to 
the AGN luminosity) is taken to be the inflow rate into a radius 
< 0. 1 pc. This radius approximately represents the radius where we 
expect the traditional viscous Q-disk to begin ( [Goodman|2003[ l, so 
the further inflow rate should be relatively constant at smaller radii 
(see the discussion in [Hopkins & Quataert|20i0l l. In the interme- 
diate and galaxy-scale simulations, we do not resolve the < 0. 1 pc 
radii from which we would like to determine the BHAR (except 
in our few extremely high-res experiments). We therefore follow 
[Hopkins & Quataert[ l |2010| l and "map" between scales, by taking 
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Figure 1. Top: Galaxy merger simulation, near coalescence of the two 
galaxy nuclei (at t 1.1 Gyr). The total SFR peaks near this coalescence, 
as inflows first reach ^kpc scales. The BH accretion rate (here multiplied 
by 500 for ease of comparison) grows rapidly during this time, but suffi- 
cient gas remains to fuel BH growth until ~ 10** yr later. Bottom: Zoom-in 
simulation of the central < 10 pc region of an active system near the peak 
in BH accretion. The SFR on small scales (< 20pc) again peaks near the 
beginning of the simulation, when the gas fraction in this radius is nearly 
unity and the gas has just reached small radii. Inefficient angular momentum 
transfer in the gas, and the continued large supply until most is exhausted, 
means that the AGN accretion rate does not peak until ~Myr later. 



the smaller-scale simulation most closely matched to the inner con- 
ditions in our large-scale simulation at each time (within some tol- 
erance) and assuming this is a crude proxy for the average behavior 
on small scales given the large-scale conditions (see Hopkins & 
|Quataert|2010| § 6 for details). We do not intend to represent this 
as a literal "zoom-in" or exactly equivalent to using a single very 
high-res simulation. But using only our ultra-high resolution simu- 
lations (which obviate the need for this mapping) does not change 
our key qualitative conclusions (but is much noisier and covers a 
much smaller dynamic range, since there are only a few such sim- 
ulations). 

3 RESULTS 

Figure [T] shows an example of two such simulations, with the 
BHAR and SFR evaluated on different scales. First, a typical gas- 
rich major merger simulation (/gas ~ 0.4 at time of merger, between 
equal-mass otherwise Milky-Way like disks). We isolate the time 
near the final coalescence of the two galaxy nuclei, and plot the 
total SFR within the galaxy and the estimated BHAR. There is a 
clear offset. The merger drives strong gas inflows into the central 



kpc which drives the SFR to a peak at just about the time of coa- 
lescence. During this time, BH growth is accelerating rapidly, but 
there is more than enough gas to continue rapid BH growth even 
after the SFR turns over (i.e. the gas begins to be exhausted). The 
BHAR finally turns over about 10** yr later. Second, we show one of 
our nuclear-scale zoom-in simulations, appropriate for the central 
regions of this simulation near the peak in the BH activity (ini- 
tial /gas = 0.8 at small radii with a 3 x lO' Mq BH: NfShlcldens 
in |Hopkins & Quataert|2010| . We compare the SFR within some 
small resolved radius (< 20pc here), where the inflows have just 
reached (hence the large initial gas fraction), to the BHAR (gas 
inflow rate into the central < 0.1 pc region). The SFR peaks im- 
mediately (more appropriately it would have peaked even slightly 
before the beginning of this simulation), as it simply traces the gas 
supply. But the inflows to small radii take some time to become 
maximal, ~ 1 Myr here. 

Figure[2]extends this to our ensemble of simulations, plotting 
the delay At between the peak of AGN activity and the peak of star 
formation activity within a given radius R. The points at 0. 1 — 10 pc 
use just each nuclear-scale simulation individually; those at 100 — 
500 pc use the intermediate-scale simulations (including ultra-high 
resolution runs on this scale which do not need another "zoom in" 
iteration to resolve the BHAR); those at > 500pc use the galaxy- 
scale simulations. In each case we measure the total SFR within 
the annulus plotted, and compare it as a function of time to the 
estimated BHAR, each smoothed over a timescale corresponding 
to a dynamical time at the radius R (since both are strongly time- 
variable).^ 

There is a clear delay which is larger when the SFR is mea- 
sured at larger radii; this arises for two reasons. First, angular mo- 
mentum loss in gas from gravitational instabilities occurs via trans- 
fer to the stellar material, when the magnitude of asymmetries in 
the potential is sufficient to induce shocks and dissipation in the gas 
( [Hopkins et aLp 009 ; Hopkins & Quataert 2010 201 1 ). In the limit 
of a pure gas system, there is no coUisionless component to ab- 
sorb the angular momentum, and so the leading-order angular mo- 
mentum loss is reduced to second-order resonance effects ^Kalnajs| 
|I97I[ l or an effective (turbulent) viscosity f Gammie|2001^ . To lead- 
ing order, the efficiency of angular momentum exchange in some 
annulus scales as oc ( 1 — ./gas ) ; where /gas = A/gas / (-^eas + Mt) is the 



Inflows 



local gas fraction in the disk i Hopkins & Quataert|201 1 1 
that "pile up" high /gas in ~ fdyn in simulations therefore tend to stall 
until /gas can be lowered (the longer gas exhaustion timescale). 

Second, it requires much less gas to fuel an AGN, as opposed 
to a starburst. The SFR inside an annulus R peaks when the gas 
mass inside that radius is maximum. In contrast, the total mass 
needed (in principle) to fuel the BH is very small. Even if "efficien- 
cies" of gas getting to the BH are low, once sufficient fuel exists to 
grow the BH at the Eddington limit, having more gas at large radii 
will not further increase the BHAR. For a typical starburst, the gas 



^ The points which correspond to the outermost scales in each tier of 
simulations should be regarded with some caution, as the characteristic 
timescales shown could be influenced by how the simulations were re- 
scaled (since the outer boundary conditions are not necessarily identical to 
a single, arbitrarily high-resolution simulation). However, the inner regions 
can be reliably evolved for many dynamical times. We have also compared 
the values plotted with the results from our more limited suite of extremely 
high-resolution intermediate-scale simulations that "bridge" multiple scales 
and obviate the need for rescalings; these are consistent in all cases. 
^ Spherical distributions of material such as the bulge and halo enter in 
separable fashion. 
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mass within ~ 500pc at tlie time of peak SFR is > 100 times tliat 
needed to fuel the BH. It is therefore at least possible that the BH 
could continue to grow at or near its Eddington limit (so the lumi- 
nosity will continue to increase exponentially) until > 99% of the 
gas is exhausted, which for typical efficiencies could require ^100 
dynamical times. Of course, accretion is not so efficient. But in 
previous simulations of galaxy mergers with simple prescriptions 
for BH accretion proportional to either the Bondi rate ( ^Springel| 
|et al.|2005} , viscous accretion rate ( Debuhr et al.|2010| or g ravi- 
tational instability accretion rate (IHopki ns & Quataert|201 if esti- 
mated on ~ lOOpc scales, the gas supply on this scale in the late 
merger stages is so large that the estimated BHAR is at or above 
Eddington throughout the entire nuclear starburst (i.e. luminosity 
rises rapidly), until > 10** yr after the SFR peaks and begins to de- 
cline from gas exhaustion. [Hopkins & Quataert] p011| l argue that 
in the case of inflows in gas-l-stellar disk, the BHAR scales (on av- 
erage) oc 1 / ( 1 + /o //gas ) where /gas is measured on large scales and 
/o ~ 0. 1 depends on the details of the system - the important point 
is that this implies accretion relatively independent of /ga.s in the 
large-/gas limit, but then cuts off when /gas 1 . 

Both these arguments imply that the critical timescale that de- 
termines the delay between peak SFR and AGN activity is the gas 
exhaustion time. The Kennicutt-Schmidt relation implies an ex- 
haustion time Af(i?) 



(1) 



where E*, Sgas, ?dyii = i /Q, = R/Vc are the SFR density, gas sur- 
face density, and orbital frequency, and e is the efficiency (fraction 
of gas turned into stars per dynamical time). Observations suggest 
values of e ~ 0.006 — 0. 1, approximately independent of local den- 
sity ( [Krumholz & Tan|2007[[Bigiel et alj2008i|Leroy et al.|2008f . 





sured Af(i?) with a crude estimator of the dynamical time, 
0.5Myr(^/pc). This is the expectation from our above deriva- 
tion if the system had a flat rotation curve with Vc ~ 100 kms^', 
and e ~ 0.02. Specific values will depend on the galaxy prop- 
erties as well as the adopted star formation law (related to the 
uncertain sub-grid SF law). More directly, we therefore com- 
pare with the gas consumption time f(Gas Consumption), where 
?(Gas Consumption) is defined as the time from when /gas peaks 
within a given annulus (about the peak SFR time) until it falls to 
/gas < 0.2 (the exact value is arbitrary so long as /gas ^ 1, we 
choose this because it gives about the right normalization). We 
measure A//r(Gas Consumption) for each simulation, and plot the 
distribution of this ratio as a function of R. There is (unsurprisingly) 
large scatter at larger radii especially, but this provides a decent ap- 
proximation to the simulation results. 

4 DISCUSSION 

Simulations of AGN fueling by gravitational instabilities naturally 
produce a delay between the time when SFRs peak inside of a 
given annulus and the time when AGN activity peaks. This offset 
scales as the gas consumption time, ~ 10 — 100 dynamical times. 
On small scales (< lOpc), this is characteristically ~ lO'yr, ris- 
ing to a few ^ 10** yr on kpc scales. These offsets are similar to 
the magnitude of time offsets suggested by various observations on 
both small and large scales ( [Davies et al.|2007{ [Schawinski et al.| 
[2009l|Wildetal.|2010| l. 

Of course, both the AGN luminosity and SFR are highly time- 
variable (see e.g. jLevine et al.|2010f . Therefore this should be taken 
as a statement only about time-averaged accretion rates and SFRs 
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Figure 2. Top: Delay time Af, between the peak of gas inflow and star for- 
mation into a given scale (< R), and the peak of AGN activity (inferred 
inflow rate into < 0.1 pc), measured in a suite of hydrodynamic simulations 
with resolution ranging from 0. 1 pc (leftmost two points) to 3 pc (middle 
two) to 50pc (right two). The blue line compares a timescale crudely sim- 
ilar to ' f[]yn where e 1/50 (the SF efficiency per dynamical time) and 
'dyii ~ R/100kms~' are assumed. Filled circles show the median Af in a 
sample of simulations, open circles correspond to the scatter (the 25 — 75% 
interval). Bottom: Same, but for the ratio A;/f(Gas Consumption), where 
((Gas Consumption) is defined as the time between when the gas fraction 
inside R peaks (usually here at /gas = 0.8 — 1) from inflows from larger 
scales or initial conditions, and when it falls below a value /gas < 0.2. 



on a timescale of at least a few dynamical times. There will be 
considerable variation from one system to another. And even in 
a time-averaged sense, the simulations show a range of behaviors, 
with the magnitude of apparent "time offsets" varying from ~ 0. 1 — 
10 times the gas consumption time (on large scales, from ~Myr 
to Gyr timescales); there are even some systems where the AGN 
luminosity peaks before the SFR. 

Nevertheless, this may be sufficient to account for current ob- 
servations without invoking poorly-understood stellar fueling (or 
feedback) processes. Those processes may still be important, of 
course, and in future work we will investigate how the explicit 
stellar evolution models presented in Hopkins et al.| ( |2011a|c|bl > 
alter the dynamics discussed here. Observationally testable con- 
sequences of these models include the magnitude of time offsets, 
their scaling with radii and the dynamical time of the galaxy, and 
the scaling of AGN activity with gas fraction inside some radius. 
Since it is gas fraction that drives these effects, one does not expect 
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to see evidence for AGN - on average - reaching tiieir peak lumi- 
nosity while the gas fraction on large scales is very large /gas ~ 1, 
but rather expects to see that the "delayed" peak in AGN activity 
corresponds to systems with /gas < 0. 1 — 0.5. And there is some ob- 
servational evidence from metal enrichment and stellar age profiles 
that high gas densities and star formation proceed "outside-in" in 
similar fashion to that predicted here in the early stages of mergers 
I Soto & Martin|2010||Kewley et al.|20T0{[Rupke et al.|2010[|ToiTeyl 



let al.|2011| >. 

This also has important implications for AGN feedback mod- 
els. AGN luminosities in realistic hydrodynamic models do not 
peak at the same time as the gas density on small scales, but rather 
at somewhat later times when most of the gas has been exhausted 
by star formation and/or expelled by stellar feedback. To the extent 
that AGN feedback assists in the "blowout" of gas (shutting down 
star formation), it primarily sweeps up the "trickle" of remaining 
gas that would provide for low-level star formation over the next 
few Gyr, not the bulk of the galaxy gas supply (see e.g. |Hopkins| 
|et al.|2008c|b[ l. For the simulations shown here, it has been shown 
previously that the AGN feedback actually acts on a very small 
residual amount of gas, and has a nearly negligible effect on the 
peak-SFR nuclear starburst (see |Cox et a l. 2006; R obertson et al.| 
|2006a[ [Hopkins et al.|[2008a| l. Similar conclusions have recently 
been obtained using entirely different accretion and AGN feedback 
models f Pebuhr et al.|2010l l. But this may be needed to explain the 
rapid color evolution of galaxies through the "green valley" (see 
e.g. |Kaviraj et al.|2011[ i. 

This also means that galaxies at their peak in star formation ac- 
tivity may not necessarily be ideal places to look for the effects of 
AGN feedback. Indeed, observations of such systems have found 
mixed results (see e.g. |Veilleux et al.|[2005[ |Rupke et al. 
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Krug_et_ay2010| and referen ces therein). Post- starburst ( Tremonti] 
et al.|2007[|Wild et a l.'2009"C oil et al.|MT) or later -stage, AGN- 



dominated mergers j S turm et al.|2011 [Humphrey et a l .2010, D unn| 
|et al.[2010[|Bautista eral.[2010||Feruglio et al.[2010| l might be more 
promising. 

We thank Vivienne Wild, Kevin Schawinski, Ric Davies, and 
Tim Heckman for helpful discussions that inspired this manuscript. 
We also thank the anonymous referee for a number of thoughtful 
suggestions. 
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